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Unusual chemical bonds are proposed. Each bond is almost covalent but is characterized by the 
thread of a small radius, ~ 0.6 x 10“^^cm, between two nuclei in a molecule. The main electron 
density is concentrated outside the thread as in a covalent bond. The thread is formed by the electron 
wave function which has a tendency to be singular on it. The singularity along the thread is cut 
off by electron “vibrations” due to the interaction with zero point electromagnetic oscillations. The 
electron energy has its typical value of (1 — 10) eV. Due to the small thread radius the uncertainty 
of the electron momentum inside the thread is large resulting in a large electron kinetic energy 
~ 1 MeV. This energy is compensated by formation of a potential well due to the reduction of 
the energy of electromagnetic zero point oscillations. This is similar to formation of a negative van 
der Waals potential. Thread bonds are stable and cannot be created or destructed in chemical or 
optical processes. 

PACS numbers: 31.10.+Z, 31.15.A-, 03.65.-w 


I. INTRODUCTION 

Various types of chemical bonds in molecules are well 
studied. See for example [H-Q- This is a branch of chem¬ 
ical physics with the well established set of basic phe¬ 
nomena. Methods of quantum chemistry are effectively 
applied to develop that field. The common feature of 
chemical bonding is the typical energy. Despite a variety 
of bond types there is the upper limit of energies involved 
into the phenomena. It is the atomic scale energy of the 
order of ten electron volts. Also the typical spatial scale 
in chemical bonding is no shorter than the atomic size. 

It looks unusual if energies, involved into molecular 
bonding, are of the order of 1 MeV and the spatial scale 
is 10“^ of the size of hydrogen atom. At the first sight, 
it is impossible since chemical bonding is associated with 
the atomic processes. 

The point is that the electron state can have a ten¬ 
dency to be singular along the line (thread) connecting 
two nuclei in a molecule. The origin of such singularity 
is clear from Schrodinger equation in cylindrical coordi¬ 
nates which can have the solution logarithmically singu¬ 
lar on the axis. The electron “vibrates” due to the inter¬ 
action with zero point electromagnetic oscillations and 
the singularity smears out within the thread of a small 
radius 0.6 x 10“^^cm. This radius can be determined 
accounting for electron-photon interaction in quantum 
electrodynamics since due to short scales the electron 
momentum becomes relativistic. 

The exact state, where the electron is coupled to pho¬ 
tons, is steady and is characterized by the certain total 
energy which is conserved. In terms of Schrodinger equa¬ 
tion it would an energy eigenvalue. The state, considered 
either far from the thread or inside it, corresponds to 
the same energy. The large contribu tion to the electron 
kinetic energy y' {mc^Y + ~ — 2.93 MeV 

inside the thread should be compensated by a reduction 
of the electromagnetic energy. The latter can be roughly 


interpreted as a potential well at the thread region. One 
can say that the thread state with the chemical electron 
energy (1 — 10) eV corresponds to upper energy states in 
the deep potential well, of the order of 1 MeV, and of the 
radius 0.6 x 10“^^cm. This well is extended along the 
thread. The electron number is determined by the outer 
region with respect to the thread. This region is of the 
order of the Bohr radius as in a covalent bond. 

This is not the unique situation of application of quan¬ 
tum electrodynamics in chemistry. The famous example 
is van der Waals forces when quantum electrodynamics 
is involved. In that case the certain attraction poten¬ 
tial is also formed due to the reduction of the energy of 
electromagnetic zero point oscillations. 

The high energy parts, involved into thread formation, 
correspond to the typical time 10“^^s. Optical processes 
are slow compared to that time. They lead to an adi¬ 
abatic motion of the thread parameters but absorption 
probability of such quanta is exponentially small. The 
life time of the thread bond itself is exponentially large 
(practically infinity) until a high energy particle or 7 - 
quantum destroys it. Also the thread bond cannot be 
created in chemical or optical processes. The above high 
energy impact is necessary for that. 

Thread bonds can be created by an external high en¬ 
ergy radiation. One can put a question on a role of such 
stable bonding in biological processes. 

In section m the singular (between two Coulomb cen¬ 
ters) solution of Schrodinger equation is obtained. In 
section mil the inner thread properties are studied. 


II. MODEL WITH TWO COULOMB CENTERS 

Electrons wave function in a molecule is complicated. 
The electron motion is collective and is not reduced to 
a single electron problem. The electron wave function 
is smooth in the space. We focus on the certain speci- 
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ficity of this function. If to consider usual atomic sizes 
the electron wave function can be formally singular on 
some line connecting two nuclei in the molecule. In the 
exact wave function the singularity is cut off at distances 
shorter than atomic ones. 

Therefore the problem is separated in two parts. First 
one should find the singular solution of quantum mechan¬ 
ical equations for electrons in the molecule. Second one 
should analyze mechanisms of smearing of that singular¬ 
ity. In this section we focus on the singularity formation 
on the line connecting two nuclei which can be treated 
in the molecule as static ones according to the Born- 
Oppenheimer approximation. 

Instead of solving the problem in full one can simplify 
it considering a single electron. The main features of 
the singular wave function can be demonstrated studying 
the artificial situation when one electron is in the field of 
two positive charges fixed at two points. The singularity 
line, connecting two nuclei in the molecule, is associated 
with large energies and is hardly influenced by chemical 
forces. For this reason the singularity line between two 
real nuclei in the molecule is of the same type as between 
two point charges formally fixed in the space. 


A. Formalism 


To study the electron in the Coulomb field of two pos¬ 
itive point charges at the points z = ±(t it is convenient, 
instead of cylindrical coordinates r = ^ + 2 /^, and 
ip, to use the elliptic ones 77 , and p [J 


c.»? = 


r2 ± ri 
2a 


( 1 ) 


Here ri ,2 = -\/(z =F distances to the Coulomb 

centers shown in Fig. [T] The surface of a constant ^ is 
the ellipsoid 



U 2 (e - 1 ) 


( 2 ) 


with the focuses at z = ±cr. The surface of a constant 77 
is the hyperboloid 

^_1_ = 1 ('Q'l 

aW 0-2(1-772) 

with the focuses at the same points. The coordinate ^ 
takes values from 1 to 00 and 77 from —1 to 1. Intersec¬ 
tions of the surfaces ([ 2 ]) and ([3|) with the plane 77 = 0 are 
shown in Fig. [TJ 

The Coulomb interaction potential is 


^ ^ 2Ze^ C 

ri r 2 cr ^2 _ jy2 ’ 


(4) 


where Ze is the positive charge at each center. The 
Schrodinger equation for the electron 
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- + UiI; = EiIj 

2m 


( 5 ) 



FIG. 1: Intersections of surfaces of constant elliptic coordi¬ 
nates 5 and 77 with the plane y ~ 0 are shown. Two Coulomb 
centers are at the points z = icr. 


in elliptic coordinates takes the form 


^2 


27770-2(^2 _ ^2^ 

2Ze2 C 


9 f.2 , 9 2 d'^ 


a — 77 " 


Tp = Eip. 


( 6 ) 


We consider an axially symmetric wave function. Since 
the variables are separated it has the form 7 /;(^, 77 ) = 
t/’i (C) V ^2 ( 77 ) ■ One can introduce dimensionless parameters 
e = —2ma‘^E/h^ and p = AZajrB where rs = h^/{me^) 
is the Bohr radius. 

After separation of variables the two Schrodinger equa¬ 
tions are 


~ 9 - (7) 

= ( 8 ) 

where /3 is some constant. 

The variable 77 can be written as 77 = cosy. Then 
Eq. ([ 8 ]) has the form 

(siny^^ = (/3 -£:cos2x)7/72. (9) 

siny dx\ ox J 

The solution, non-singular at y = 0, should be also non¬ 
singular after continuation to y = tt. This is the condi¬ 
tion to choose the parameter /3 for a given e. When two 
Coulomb centers coincide (tr = 0 and therefore e = 0) 
that condition turns to /3 = —l{l + 1) as for Legendre 
polynomials In that case the variable y coincides 
with the azimuthal angle 9. 


B. Close Coulomb centers 

Below two close Coulomb centers are considered un¬ 
der the condition a tb- We study the state which is 
isotropic {I = 0) in the limit cr = 0. At a finite a there is 
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the small correction to the wave function which can be 
written as '02 = 1 + <5'02- Using the relation 

sinx^^ =/rfx(/3 - ecos^x)sinx (10) 

one can obtain 

S'4’2 = (I “ (2cos^ 2 ) ~ I 

The solution m is finite at x = 0 (?? = 1). In order to 
get it finite at x = t’" (^ = ~1) it should he 13 = e/3. 
Note that e ^ p'^ ^ /r^ are small. 

At tr <C rs there is a small region of the size a around 
the centers in Fig. [TJ This region corresponds to ^ ^ p ^ 
1. The region of the order of the Bohr radius is much 
larger and relates to large We consider first this region. 
In Eq. 0 one can omit /3 and to write — 1 ~ In 
the limit 1 ^ Eq. o takes the form 


52-01 2 501 p 


-£■ 01 , 


( 12 ) 


which coincides with the radial Schrodinger equation 
with Z = 0 in the Coulomb field of the point charge 
2Ze 0 . The solution of dED, decaying on infinity and 
finite at small distances, corresponds to the eigenvalue 
£ = p2/4. This value relates to the ground state energy 
in the Coulomb field of the point charge 2Ze. 

For our purposes one needs a solution which also de¬ 
cays on infinity but is singular at r = 0 and <0-2. To 
obtain that one should write the solution of (1121) in the 
form 0 

V'l(C) = C^U 2 ^/^)-l 

where 


G{a, I3,v) — l-\- — -I- 


a/3 a(a + l)/3(/3 + 1) 


2 !i ;2 


(14) 


To obtain the solution at shorter distances one should 
use the asymptotics of the function m at n <C 1. Thus 
at small [1 —p/{2y/e)] one can obtain from Eq. (fT^ @ 


27 ?) 1 «««’■«/<’■ (15) 

On the other hand, at not large ^ the left-hand side of 
Eq. ([71) is the principal one and the solution is 

( 16 ) 

Eq. (fT51) goes over into the form (|T5]) when their applica¬ 
bility intervals overlap. 


The wave function along the line, connecting two 
Coulomb centers in Fig. [TJ now can be written at — 
cr^ I, <C in the form 


0 (r, z) = 1 - 


rs 

16Za 



X In 


_ 8 u 2 _ 

z^ — cr^ -I- (z^ — cr^)^ -I- dcr^r^ ’ 


(17) 


where Eg = —m{2Ze^Y/{2h'^). Eq. (fT71) is also valid in 
the vicinity (r^ ^ (cr^ — z^)) of the entire line (z^ < ct^) 
between the centers, where 


'0(r, z) = 1 - 


tb 

8Za 


E \ , 2Vcr2 - z2 


1 -In 

Eo 


(18) 


At large distances, as follows from (HI - da and the 
definition o, 


Vi(r, z) = rV^o-1 exp , rs « R, 

(19) 

where + z^. 

The wave function of the electron in the Coulomb field 
of two positive point charges Ze exponentially decays at 
large distances (HU). But on the line, connecting two 
charges, 0 has the logarithmic singularity (fT51) if the en¬ 
ergy does not coincide with the eigenvalue Eq. The ab¬ 
sence of singularities is a usual condition to determine an 
eigenvalue. The eigenvalue Eq coincides with one in the 
Coulomb field of one point charge 2Ze. In the limit con¬ 
sidered, a <C rs, corrections to that eigenvalue are small. 
Eq is the lowest energy level shown in Eig.|2](b) where the 
potential energy (|4|) relates to the curve with the dashed 
insert. The eigenfunction is plotted in Fig. UKa) with the 
dashed part of the curve. 

When the energy differs from Eq the singular wave 
function is not physical at the first sight. But the situa¬ 
tion is more complicated as described in Sec. IIIIl 


III. ELECTRON THREADS 


The electron wave function, for the artificial case of 
two static point charges, and the wave function for two 
partner nuclei in a real molecule have common proper¬ 
ties. This happens since at short distances (less than tb) 
atomic forces are not significant. In the real molecule 
close to a nucleus, that is below the shell of inner elec¬ 
trons, the wave function is of the type H 3 - The same is 
valid for the nucleus at a partner site. Between the part¬ 
ner sites the wave function is logarithmically singular, as 
(HS, along the connecting line —<7<z<a. In a real 
molecule cr ^ r^. In a close vicinity (r <C cr) of the line, 
connecting two partner sites, the wave function, with the 
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FIG. 2: Features of the thread state at j/ = 0 and constant 
r}. Dashed curves correspond to the absence of tread. Thin 
vertical lines restrict the thread region, (a) Thread wave func¬ 
tion. (b) Energy levels in the potential well in the absence of 
thread. 


logarithmic accuracy, can be written in the form 

f A{a) ln [-\/{z — aY +r'^ + {z — cr)], \z — a\ a 
'i/’ = ■! 2A(z) In r , \/— z^ ^ a 

[ A(—cr) ln [-\/(z + cr)2 + — [z + a)], |z -|- cr| <C cr 

( 20 ) 

Here A{z) is the certain function accounting for real in¬ 
tramolecular forces. The exact form of this function is 
not crucial for our purposes. 

An electron, acted by an external macroscopic field, is 
spatially smeared due to the interaction with zero point 
electromagnetic oscillations @. If ft is the electron dis¬ 
placement from its mean field position then (u) = 0 and 
i^T = \/{u^) is finite. In an external field (for example 
Coulomb one) due to the uncertainty in electron positions 
it probes various parts of the potential. This changes its 
energy. The famous example of this phenomenon is the 
Lamb shift of atomic levels 0. 

The electron momentum is not conserved in virtual 
processes and the electron “vibrates”. The uncertainty 
tt results in small corrections (such as Lamb shift) when 
the spatial scale of the electron wave function is larger 
than Tt- But when the wave function, as a solution of 
quantum mechanical equations, has a singularity that un¬ 
certainty plays a significant role resulting in smearing of 
the singularity within ry. 

Therefore the singularity line between two nuclei in 
a molecule has the finite width of the order of r^. It 
can be called thread bond. In Fig. [2](a) the wave func¬ 
tion increases towards the tread and it is cut off on the 
thread size r^. The potential energy in Fig.j^Kb) is of the 
usual Coulomb type outside the thread region. Inside the 
thread the mean field potential no more exists because 
photon degrees of freedom enter the game. 

The singularity of the wave function is cut off on 
the distance tt and the state becomes physical. So 


the expression (l20t corresponds to not small distances 
r^, |(j^ — z'^\. According to Eqs. (IT71) and (ITOl) . the sin¬ 
gular solution of the Schrodinger equation exists at any 
pre-logarithmic factor (1 — E/Eq). This means that the 
electron energy of the thread state is not quantized due 
to the interaction with photons. This can be explained 
in another way. Broadening of discrete levels in a well 
(quasi-continuous density of states) is narrow and can be 
calculated from the perturbation theory. Calculation of 
smearing of the singular wave function requires all or¬ 
ders of the perturbation theory. In this case the density 
of states is smoothly continuous. 


A. Estimate of the thread radius 

Unlike the Lamb shift, it is impossible to apply the 
perturbation theory to calculate smearing of a singular 
wave function of electron. In that problem one should 
account for all orders of the perturbation theory. For 
this reason, we use below the approximate method just 
to estimate the electron mean squared displacement Q. 
The method is successfully applied for study of the Lamb 
shift. In this method the electron motion under the ac¬ 
tion of zero point oscillations can be described by the 
equation 


d^u 

m—TT + = —e£, (21) 

dE 

where m is the electron mass and U ^ me'^/h^ is the 
electron rotation frequency in the atom. 

One can use the Fourier expansion 

u{R, t) = Uk exp{ikR — iujkt) (22) 

k 

and analogous one for the fluctuating electric field 
£{R,t). Since u{R,t) is real it should be = U-k and 
LO-k = —tOk in the expansion The condition uk 
has to be held in this method. The solution of Eq. (ED 
is of the form 

1 + 1 

2Tfi\ujk I \ \a)k I “t" H IWfcI U 

According to the quantum mechanical approach, Eq. (1231) 
should be modified as 

^ _ e£k f yi + rik \fnk \ 

2TO|wfc| V |wfc| + H |wfc|-U/’ 

where rik is the number of quanta, the first term relates 
to the quanta emission, and the second one to the ab¬ 
sorption. 

The mean squared displacement is 

{n^)= f ( 25 ) 

k 
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where V is the volume of the system. Since in our case 
rife = 0 , the mean squared displacement has the form 


(u^} = 




W) 


Am? ^ w^dwfel +f2)2 


(26) 


The energy of zero point oscillations is 


At: 


It follows from here that (|ffep) = 27r?i|a;fc|/T^. Using the 
summation rule 



Airk'^dkV 

(27r)3 


(28) 


(the coefficient 2 accounts for two photon polarizations) 
and the relation uik = ck, one can obtain from Eq. (ESI) 

a 



TT 



ujdu! 

(w +12)2 


(29) 


The upper limit ujmax is determined by the condition of 
non-relativistic motion, that is ujmax — m?/h. In the 
relativistic region Uk decreases due to the enhancement 
of the relativistic mass. Using the above rydberg estimate 
for hn one can obtain for the thread radius 


rr = V (m^) = In ^ > (30) 

where Tc = h/{mc) ~ 3.86 x 10“^^cm is the electron 
Compton length. The exact cut off 12 is not crucial for 
the estimate. The applicability condition of the approach 
used (u^)fcmaa; “C 1 is valid since (u^)/r^ ~ 

e^/ {he). 


B. Thread properties 

When the wave function is smooth one can consider 
the electron motion in the effective potential {V{R + u)) 
averaged on fast electron motions u. Expansion up to the 
second order in u and using Eq. (l30l) produces the Lamb 
shift 0 coincided (excepting the numerical coefficient) 
with the exact result 

When the solution of quantum mechanical equations 
is singular, it is cut off at the distance rr- The thread 
radius can be estimated from Eq. (1501) as rr — 0.15 rc — 
0.58 X 10“^^ cm. 

At distances tt < r radiation effects are small and 
one can use the quantum mechanical description of the 
electron. At rc < r this the Schrodinger formalism but 
at rx < r < Tc one should apply Dirac equations. In 
this case, besides the term Inr, the wave function also 
contains the term r^jr which is small at rc <C r. See 


Appendix (|A.4I) . The total number of electrons is deter¬ 
mined by the region r ^ tb outside the thread as in a 
covalent bond. The fraction of electrons in the region 
r < rc is small as r^/r^ ^ {e^ jhc)^. 

The wave function is schematically plotted in Fig.j^Ka) 
where it almost corresponds to the quantum mechanical 
approach (weak radiation corrections) outside the thread. 
Inside the thread the wave function has the peak which 
can be interpreted as one resulted from a narrow poten¬ 
tial well. Inside the thread the electron state is coupled to 
the electromagnetic subsystem and cannot be considered 
separately. Only the total thread energy has a meaning 
without its separation on electron and electromagnetic 
ones. 

The exact state, where the electron is coupled to pho¬ 
tons, is steady and is characterized by the certain total 
energy which is conserved. In terms of Schrodinger equa¬ 
tion it would an energy eigenvalue. The state, considered 
either far from the thread or inside it, corresponds to the 
same energy. The large contribution to the electron ki¬ 
netic energy + {hcjrT)'^ — m? ~ 2.93 MeV in¬ 

side the thread should be compensated by a reduction of 
the electromagnetic energy. The latter can be interpreted 
as a potential well at the thread region in Fig. [IKb) in¬ 
stead of the dashed curve indicating the usual Coulomb 
potential. One can say that the thread state with the 
chemical energy scale E ^ {\ — 10) eV corresponds to 
upper energy states in the deep potential well, of the or¬ 
der of 1 MeV, and of the radius rr- This well is extended 
along the thread. 

We emphasize that the interpretation in terms of the 
potential well is approximate since this is not a single 
particle quantum mechanics but coupling to the electro¬ 
magnetic system. 

The reduction of the electromagnetic energy reminds 
the van der Waals phenomenon when the negative poten¬ 
tial between atoms is also formed by the reduction of the 
energy of zero point electromagnetic oscillations ii- 

The electron thread energy is continuous and when it 
is smaller than the lowest level in Fig. IHb) the thread 
state can be destroyed by quanta absorption only. The 
high energy parts, involved into thread formation, cor¬ 
respond to the typical time I 0 “^ 2 g Optical processes 
are slow compared to that time. They lead to an adia¬ 
batic motion of the thread parameters but the absorption 
probability of such quanta is exponentially small Q . The 
life time of the thread bond itself is exponentially large 
(practically infinity) until a high energy particle or 7 - 
quantum destroys it. Also the thread bond cannot be 
created in chemical or optical processes. The above high 
energy impact is necessary for that. 

The length of the thread bond and a form of the outer 
electron cloud are determined by chemical mechanisms 
of interaction with surrounding electrons. The electron 
density is mainly localized outside the thread on the usual 
distance of the order of the Bohr radius as in a covalent 
bond. 
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FIG. 3: Schematic representation of molecular covalent 
bonds. The thread bond between two nuclei is shown by the 
line. 

IV. DISCUSSIONS 

It looks unusual if energies, involved into molecular 
bonding, are of the order of 1 MeV and the spatial scale 
is 10“^ of the size of hydrogen atom. At the first sight, 
it is impossible since chemical bonding is associated with 
the atomic processes. The point is that the electron state 
can have a tendency to be singular along the line (thread) 
connected two nuclei in a molecule. This singularity is 
cut off by electron interaction with zero point electromag¬ 
netic oscillations. The electron also “vibrates” within the 
thread of the small radius tt where the singularity smears 
out. 

The singularity smearing is associated with non¬ 
conservation of the electron momentum. The uncertainty 
of this momentum is h/ y/JvF). In the absence of exter¬ 
nal potential, as follows from the mean squared dis¬ 
placement is infinite (17 = 0). In this case the uncertainty 
of the electron momentum is zero and the singularity is 
not cut off. This also follows from the exact spectral 
representation of the electron Green’s function in the ab¬ 
sence of external fields Q. 

The peak of electron wave function in the thread re¬ 
gion is associated with formation of the certain deep 
1 MeV) potential well along the thread. The po¬ 
tential well is due to the reduction of the energy of elec¬ 
tromagnetic zero point oscillations. The interpretation 
of thread bonding can be in terms of an electron state 
with the energy E {1 — 10) eV which corresponds to 
upper energy states in the deep well. The interpretation 
in terms of the potential well is approximate since this is 
not a single particle quantum mechanics but coupling to 
the electromagnetic system. 

This phenomenon is described by quantum electrody¬ 
namics since due to short scales the electron momentum 
becomes relativistic. This is not the unique situation 
of application of quantum electrodynamics in chemistry. 
The famous example is van der Waals forces when quan¬ 
tum electrodynamics is involved. In that case also the 
certain steady attraction potential is formed due to the 
reduction of the energy of electromagnetic zero point os¬ 


cillations. 

Since high energy is involved into the thread formation 
it cannot be created (or destroyed if it exists) by chemical 
or optical processes. It can be created by a high energy 
particle or y-quantum. After that the thread bond lives 
the exponentially large (practically infinity) time. A high 
energy process, of the above type only, can destroy it. 

The schematic representation of the thread bond in 
a molecule is in Fig. [31 The thin thread connects two 
partner sides. The thread bond is associated with the 
outer covalent part where the electron cloud is mainly 
localized (not on the thread). 

External high energy radiation (natural or artificial) 
can result in accumulation of thread bonds in biological 
molecules, for example in DNA. One can put a question 
how thread bonds will affect biological processes, for ex¬ 
ample, replication of DNA. 


V. CONCLUSIONS 

Unusual chemical bonds are proposed. Each bond is 
almost covalent but is characterized by the thread of a 
small radius, ~ 0.6 x 10“^^cm, between two nuclei in 
a molecule. The main electron density is concentrated 
outside the thread as in a covalent bond. The thread is 
formed by the electron wave function which has a ten¬ 
dency to be singular on it. The singularity along the 
thread is cut off by electron “vibrations” due to the inter¬ 
action with zero point electromagnetic oscillations. The 
electron energy has its typical value of (1 — 10) eV. Due 
to the small thread radius the uncertainty of the elec¬ 
tron momentum inside the thread is large resulting in a 
large electron kinetic energy ^ 1 MeV. This energy is 
compensated by formation of a potential well due to the 
reduction of the energy of electromagnetic zero point os¬ 
cillations. This is similar to formation of a negative van 
der Waals potential. Thread bonds are stable and cannot 
be created or destructed in chemical or optical processes. 
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Appendix: SINGULAR SOLUTION OF DIRAC 
EQUATIONS 

Schrddinger equation formally has the solution Inr in 
cylindrical coordinates. Below we establish the continu¬ 
ation of this singular solution to the region r < where 
one should use the Dirac formalism. In this case the 
wave function is bispinor consisting of two spinors (p and 
X @. Since we are interested by the singular wave func¬ 
tions (large kinetic energy part) one can ignore potential 
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energy and consider free electron Dirac equations 

[Eq — ihc{aS/)\ X = mc^(p , [Eq + jSc((tV)] tp = rnc^x^ 

(A.l) 

where Eq is the total relativistic energy and a are the 
Pauli matrices Substituting the spinor x from the 
second equation (lA.ll) into the first one we obtain 

[Eq + (fi = m^c^ip. (A.2) 

The equation for y is the same. We use here the rela¬ 
tion (tTV)(CTV) = V^. These equations have solutions 
N()(r\/En/c^ ^ m^c^/h) and {aV)No where Nq is the 
Neuman function [lOj . 

We use the asymptotic forms Nq ^ \nr and VfVo ~ 
r/r^. Those solutions are compatible with the equations 


(lA.ll) if to put 

ZT ZT 

if = folnr- ■^{ra)ifo, X = V^olnr-f ■^{ra)fo, 

(A.3) 

where fo is some constant spinor. Eqs. (IA.3I) account 
for the condition that for our non-relativistic energies 
{Eq — mc^) <C mc^, two spinors should coincide at large 
distances rc r Q . In the standard representation $ = 
f + X 0 = 7^ ~ X 

<^ = 2fQhir, Q = -^{fa)fo (A.4) 

at large distances (non-relativistic limit) the wave func¬ 
tion is the usual spinor $. 
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